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Generalizing the mode-coupling theory for ideal liquid-glass 
transitions, equations of motion are derived for the correlation 
functions describing the glassy dynamics of a diatomic probe 
molecule immersed in a simple glass-forming system. The 
molecule is described in the interaction-site representation 
and the equations are solved for a dumbbell molecule consist- 
ing of two fused hard spheres in a hard-sphere system. The re- 
sults for the molecule's arrested position in the glass state and 
the reorientational correlators for angular-momentum index 
£ — 1 and £ — 2 near the glass transition are compared with 
those obtained previously within a theory based on a tensor- 
density description of the molecule in order to demonstrate 
that the two approaches yield equivalent results. For strongly 
hindered reorientational motion, the dipole-relaxation spectra 
for the a-process can be mapped on the dielectric-loss spectra 
of glycerol if a rescaling is performed according to a sugges- 
tion by Dixon et al. [Phys. Rev. Lett. 65, 1108 (1990)]. It is 
demonstrated that the glassy dynamics is independent of the 
molecule's inertia parameters. 

PACS numbers: 64.70.Pf, 61.25.Em, 61.20.Lc 



I. INTRODUCTION 

The mode-coupling theory (MCT) for the evolution of 
structural relaxation in glass-forming liquids was origi- 
nally developed for atomic systems and for mixtures of 
atoms or ions. Detailed tests of the theory have been 
provided by comparisons of the predictions for the hard- 
sphere system (HSS) with dynamic-light-scattering data 
for hard-sphere colloids, as can be inferred from Ref. J|] 
and the papers quoted there. Quantitative tests have also 
been made by comparing molecular-dynamics-simulation 
data for a binary mixture with the MCT results for the 
model • A series of general implications of the MCT 
equations were derived like scaling laws and relations be- 
tween anomalous exponents describing power-law spec- 
tra and relaxation-time scales which establish some uni- 
versal features of the dynamics It was conjectured 
that these results also apply to molecular liquids. In- 
deed, there is a large body of literature, which is re- 
viewed in Ref. || , dealing with the analysis of data from 
experiments or from molecular-dynamics simulations for 
complicated systems in terms of the universal MCT for- 
mulas. These studies suggest that MCT describes some 
essential features of the glassy dynamics for molecular liq- 
uids. Therefore, it seems desirable to develop a detailed 



microscopic theory also for systems of non-spherical con- 
stituents. 

A mode-coupling theory for molecular systems has 
been studied in Refs. [[7j-jl2fl where the structure is de- 
scribed by tensor-density fluctuations. The basic con- 
cepts of the MCT for simple systems like density cor- 
relators and relaxation kernels have been generalized to 
infinite matrices. The equations for the non-ergodicity 
parameters and critical amplitudes were solved. These 
quantities generalize the Debye- Waller factors of the ar- 
rested glass structure and characterize its changes with 
temperature. Comparison of the theoretical findings with 
molecular-dynamics-simulation data for water PJll[| and 
for a system of linear molecules |n| demonstrates that 
the theory can cope with microscopic details. However, 
the derived equations are so involved that further simpli- 
fications would be required before correlators or spectra 
could actually be calculated. 

The simplest question of glassy dynamics of the rota- 
tional degrees of freedom concerns the motion of a single 
linear molecule in a simple liquid. This problem is equiva- 
lent to the study of a dilute solution of linear molecules in 
an atomic liquid as solvent. For this system, a MCT has 
been developed, generalizing the equation for a tagged 



particle in a simple lie 
for a tagged molecule 



uid to an infinite-matrix equation 
13f| . The equations were solved for 
a molecule consisting of two fused hard spheres immersed 
in a HSS |l^,|lj]. The validity of the universal laws for 
the reorientational dynamics was demonstrated. Char- 
acteristic differences for the a-process of the relaxation 
for angular-momentum index i = 1 and i = 2 were iden- 
tified which explain the differences of spectra measured 
for I = 1 by dielectric-loss spectroscopy and for £ — 2 
by depolarized-light-scattering spectroscopy. The exper- 
imentally established large ratio of the a-relaxation times 
for the £ = 2-reorientational process and the longitudi- 
nal elastic modulus was also obtained |l4j| . These exam- 
ples show that MCT can provide general insight into the 
glassy dynamics of rotational degrees of freedom which 
goes beyond the contents of the universal formulas. 

Within the basic version of MCT, the tagged-particle- 
density-fluctuation correlator for wave number q consid- 
ered as function of time t, cf> q (t), or the dynamical struc- 
ture factor for frequency to, S q (u>), can be written as: 
<f>q(t) = <t> S q*(t/to) and S q (uj) = S*(u;t ). Here the func- 
tions 4> s q * (t) and S* (u>) are completely determined by the 
equilibrium structure. This holds for times outside the 
transient regime, t/to 3> 1, or for frequencies below the 
band of microscopic excitations, uto <C 1. The subtleties 
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of the transient dynamics like the dependence of oscil- 
lation frequencies on mass ratios enter the long-time dy- 
namics and the low-frequency spectra via a common time 
scale to only. This means that the statistical information 
on the long-time dynamics is determined up to a scale 
to by the statistics of the system's orbits in configura- 
tion space rather than by the orbits in phase space. The 
glassy dynamics as described by functions like 4> s q (t) or 
S q (uj) deals with the probabilities of paths through the 
high-dimensional potential-energy landscape. The com- 
plicated dynamics on microscopic time scales is irrelevant 
in the long-time regime; it merely determines the scale 
to for the exploration of the configuration space. The 
cited results of the MCT for simple systems and mix- 
tures [[s| |l5|]l6fl are not valid for the mentioned theories 
for molecular systems J?] |ll|,^,^4| , which imply isotope 
effects for the glassy dynamics. A change of the mass 
ratio of the molecule's constituents shifts the center of 
gravity, and thereby the mode-coupling coefficients are 
altered. This leads to shifts of the glass-transition tem- 
perature, the particle's localization lengths, and the like. 
In this respect a system of A-B molecules would behave 
qualitatively different than an A-B mixture. There are 
no experimental observations demanding that the long- 
time dynamics is independent of the inertia parameters 
of the molecules. But we consider the specified iso- 
tope effects as artifacts of the approximations underly- 
ing the so far studied extensions of MCT. This critique 
and the formidable complexity of the theories based on 
the tensor-density descriptions appear as a motivation to 
search for an alternative approach towards the glassy dy- 
namics of molecular systems. An alternative MCT was 
proposed by Kawasaki jl7[ |. But so far, nothing is known 
about the solutions of his equations nor on the results 
concerning the inertia-parameter issue. In this paper the 
suggestion of Chong and Hirata [Q will be followed, and 
the MCT will be based on the interaction-site represen- 
tation of the system [ p~9|j20| . 

The description of a molecular liquid by interaction- 
site densities is inferior to the one by tensor densities. 
The correlators of tensor densities can be used to express 
the ones of interaction-site densities but not vice versa. 
Interaction-site theories have also difficulties to handle 
reorientational correlators. Therefore, it is a major issue 
of this paper to show, that the indicated ad hoc objec- 
tions against a MCT based on an interaction-site repre- 
sentation do not fully apply if the theory is restricted to 
a parameter regime where the cage effect is the domi- 
nant mechanism for the dynamics. To proceed, the same 
dumbbell-molecule problem shall be studied, which was 
analyzed previously Jl3| , pT| . 

The paper is organized as follows. The basic equations 
for the model are "introduced in Sec. || Then, the MCT 
for a diatomic molecule in a simple liquid is formulated in 
Sec. HI. The major problem is the derivation of formu- 



theory for the dumbbell in a HSS are discussed. The 
findings are summarized in Sec. [v| 



II. THE MODEL 

A system of N identical atoms distributed with density 
p at positions r K , n — 1, • • • , N, is considered as solvent. 
The structure can be described by the density fluctu- 
ations for wave vectors q: Pg- = X) K ex P(*<? ' *«)• The 
structure factor S q — (| p q \ 2 )/N provides the simplest 
information on the equilibrium distribution of these par- 
ticles. Here (• • •) denotes canonical averaging for tem- 
perature T. Because of isotropy, S q only depends on 
the wave- number q =| q |. The Ornstcin-Zernike equa- 
tion, S q = 1/[1 — pc q ], relates S q to the direct cor- 
relation function c q . The structural dynamics is de- 
scribed in a statistical manner by the normalized den- 
sity correlators (f> q (t) — (pf(t)*p$)/NS q . They are real 
even functions of time t and exhibit the initial behavior: 

is 



<j> q {t) = 1 - \{n q tf + 0(\t\ 3 ). He re Sl v = qv/y/S~ q 

the bare phonon dispersion; v = \J ksT/m denotes the 
thermal velocity of the particles with mass m |2(J . 

A rigid molecule of two atoms A and B shall be con- 
sidered as solute. Let f a , a = A or B, denote the 
position vectors of the atoms, so that L = \?a — re| 
denotes the distance between the two interaction sites. 
Vector e = (?a — ^b)/L abbreviates the axis of the 
molecule. If m a denotes the mass of atom a, the to- 
tal mass M = tua + ms and the moment of iner- 
tia I — mATTiBL 2 JM determine the thermal velocities 
Vt = y/ksT/M and vr = \JkBT jl for the molecule's 
translation and rotation, respectively. Let us introduce 
also the center-of-mass position fc = {mA'rA+'mB'rB)/M 
and the coordinates z a of the atoms along the molecule 
axis: za = L(vhb / 'M), zb — —Hjtia/M). The posi- 
tion of the molecule shall be characterized by the two 
interaction-site-density fluctuations 



pi = cxp(iq- f a ), a = A or B. 



(1) 



The two-by-two matrix w g of static fluctuation corrcla- 

4*4) 



tions w q b = (p^pt) is given by 



w a q b =S ab + (l-S ab )j Q (qL), (2) 

where here and in the following ji(x) denotes the spher- 
ical Bessel function of index £. The solute-solvent in- 
teraction is described by the pair-correlation function 
hq = {p* q P a ^ilPi which is expressed by a direct correla- 
tion function c q jl9| 



b 



(3) 



las for the mode-coupling coefficients. This will be done 
within the Mori-Fujisaka formalism, and the details are 
presented in Appendix M. In Sec. IV, the results of the 



The dynamics of the molecule shall be characterized 
by the interaction-site-density correlators 



2 



*y 6 w = (4) 

These are real even functions of time obeying F^ b (t) — 
F ba (t). They shall be combined to a two-by- two-matrix 
correlator F q (t). Its short-time expansion can be noted 



F 9 (t) 



ig 2 J 9 i 2 + 0(|t|- 



(5) 



The continuity equation reads pi — iq ' ■ j'2, where the 
current fluctuation is j% = v a pi with v a denoting the 
velocity of atom a. Therefore, one gets J° b = ((q-j~)*(q- 

Ab 



J0) /q 2 . The result splits in a translational and rotational 



part, 3 q = J q + Jf, where |21|]: 



jTab _ 2 ab 
jRab _ „,2 ^2 



(6a) 



(l-0(io(gi)+i 2 (g£))]. (6b) 

Let us note the small-g expansion of the density corre- 
lators in the form 



F« b (t) = l-±q 2 C ab (t) + 0(q% 



(7) 



The diagonal elements of the symmetric matrix C(t) are 
the mean-squared displacements 

5r 2 a (t) = ((r a (t)-r a (0)) 2 )=C aa (t), (8a) 

while the off-diagonal elements can be related to the 
dipole correlator 

C x {t) = (e(t) • e) 

= [C AB (t) - \{C AA {t) + C BB {t))\ / L 2 . (8b) 

The mean-squared displacement of the center of mass can 
be expressed as 

= [m A Sr 2 A (t)+m B Sr 2 B (t)}/M 

+ (2I/M)[C 1 (t)-l}. (8c) 

Expanding Eq. (|5|) in q yields the initial decay 

C(t) = C + 3J t 2 + O(|<| 3 ). (9) 

Here the initial value Co is due to the expansion of 
Eq. (Q) , while the prefactor of the £ 2 -term is due to the 
zero- wave- number limit of Eqs. (^|): 



o x (q) = (p A ±p B )/V2, 



x 



Px\h i — \rq -■- rq 

The transformation matrix P = P^ 1 



N or Z. 

reads 



1 

71 



i i 

i -i 



It diagonalizes the matrices w q from Eq. 
Eqs. (§): 



(11a) 



(lib) 



and J q from 



(Pw,P)*» = i*»«;,(«), w x (q) = l±j (qL), (11c) 
(P J* P)*» w x (q), (lid) 

(P jRp r = 5-vlvlL 2 [1 T (io(gi) +ia(g£))], (He) 



where x, y = or Z . Also the matrix of density correla- 
tors is diagonalized. Introducing the normalized correla- 
tors (j)g(t), one gets 



(PF q (t)PT y = S xy cj> x q (t)w x (q). 



(12a) 
(12b) 



The mean-squared displacements are equal and shall be 
denoted by Sr 2 (t) = 5r 2 A {t) = 5r 2 B (t), so that Eq. @ 
reads 5r 2 {t) = 5r 2 c {t) + (l/2)i 2 [l - d(t)]. The matrix 
C(t) is diagonalized 



(PC(t)P)™ = 25r 2 c (t)+L 2 , 
(PC(()P) ZZ = -L 2 d(t). 



(13a) 
(13b) 



In Appendix A it is shown how the correlation func- 
tions in the interaction-site representation can be ex- 
pressed in terms of the ones in the tensor-density rep- 
resentation. 



III. APPROXIMATIONS 

A. The solvent-density correlator 

The density correlator of the solvent is needed to for- 
mulate the equations for the probe molecule. This quan- 
tity is discussed comprehensively in the preceding liter- 
ature on the MCT for simple systems f22]| . Let us note 
here only those equations which have to be solved in or- 
der to obtain the input information for the calculations 
of the present paper. First, there is the exact Zwanzig- 
Mori equation of motion ]2(J relating the correlator for 
density fluctuations 4> q (t) to the correlator m q (t) for the 
force fluctuations: 



cf = l 2 (i - s ab ), js" = 4 + § 4 z a z b . (io) dU q {t) + n 2 q <t> q {t) + o 2 / dt> m q (t - o cV<MO = o 

Jo 

For symmetric molecules one gets vtla = ma = M/2, 
I = ML 2 /4, and za = —zb — L/2. In this case, 
there are only two independent density correlators, since 
F AA {t) = F BB (t). It is convenient to perform an or- 
thogonal transformation to fluctuations of total number- 
densities pn{<i) and "charge" densities pz(q)' 



(14) 

Second, there is the approximate expression for kernel 
m q (t) as mode-coupling functional 



m q {t) =T q [(j>{t)]. 



(15a) 
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The functional F q is rederived as Eq. (Bll) in Ap- 
pendix [b|. The wave numbers are discretized to M values 
with spacing h: q/h = 1/2, 3/2, • • • ,M- 1/2. Then cf>(i) 
and similar quantities are to be viewed as vectors of M 
components (f> q (t), q = 1, • • • , M, and the functional is 



— ^ Vq,kp fk fp 
kp 



(15b) 



Third, Eqs. ([14|) and ( |15a| ) imply the equation for the 
long-time limit f q — <f) q (t — > oo): 



/ ? =^[/]/{l+^[/]}- 



(16) 



For the liquid state, there is only the trivial solution 
f q = 0. The glass is characterized by a non-ergodicity 
parameter < f q < 1, which has the meaning of the 
Debye- Waller factor of the arrested structure. At the 
liquid-glass transition, the long-time limit of the correla- 
tor changes discontinuously from zero to the critical value 
f c q > o. 



The equation for the non-ergodicity parameters of the 
molec ule, f° bo ° = F ^ b (t — > oo), can be obtained from 



Eqs. (|17a| ) and (18a). It is a matrix generalization of 
Eq. (|l£ 



^[F-ZKl + ^F 00 ,/]}- 1 - 



(19) 



If the solvent is a liquid, i.e., if f q = 0, one gets F^° = 0. 
If the solvent is a glass, the long-time limits F q ho ° can 
be non trivial. In this case, the solvent properties enter 
via the Debye- Waller factors f q , whic h renormalize the 



coupling coefficients V q cb kp in Eq. (18b) 



Let us specialize to symmetric molecules. Multiplying 
Eqs. (|l7|) to ( |l9| ) from left and right with P given by 
Eq. ( |llb| ) and inserting 1 = PP between every pair of 
matrices, all equations are transformed to diagonal ones. 
Thus, there are two equations of motion 



WW 



n 



x2 



dt'm x q (t-t')d t ^ x (t')=0, (20) 



B. The solute-interaction-site-density correlators 

For matrices of correlation functions as defined in 
Eq. (|[) , the Zwanzig-Mori formalism also leads to an ex- 
act equation of motion poj : 

d 2 t ¥ q {t) + n 2 q F q (t) + n 2 q [ dt' m q (t ~ t') d t ,F q (t') = 0. 

Jo 

(17a) 

From the short-time expansion together with Eq. (||) , one 
gets 



q 2 3 q w q 



The r.h.s. of this equation is a product of two symmetric 
positive definite matrices. Hence it can be written as 
square of a matrix tt q . Splitting off this matrix in front 
of the convolution integral is done for later convenience. 

The difficult problem is the derivation of an approxi- 
mation for the matrix m q (t) of fluctuating-force correla- 
tions, so th at the cage effect is treated reasonably. The 
result, Eq. (B17) from Appendix |^, can be formulated 
as mode-coupling functional ^F q : 



mf{t)=T?[¥{t)A{t)\ 



(18a) 



After the discretization of the wave numbers as explained 



above, J- q reads 



kp 



(18b) 



The preceding equations are matrix generalizations of the 
MCT equations for the tagged-particle-density correlator 
4> q (t) in a simple liquid [g3| . 



for x — N or Z. The two characteristic frequencies 
il q , which specify the initial decay of the correlators by 
</>-(*) = l-i(^t) 2 + 0(|i| 3 ), read: 



^ 2 = (v T qf 



Vr = (v T q) 2 



l(v R L q ) 2 [l- 3o (qL)-j 2 (qL)} 

/[1+jo («£)], (21a) 
\vl[\+j Q {qL)+ 32 (qL)]{qL) 2 

/[l- jo (?£)]• (21b) 



The relaxation kernels c an b e written as m q (t) 
^[cj) x (t), <f>(t)], where Eq. ( |B17| ) gives 



(17b) ^[f x J] = [w x (q)/q- 



dk 



(q-p/qf 



2(2tt) 3 
xw x (k)pS p c N (p) 2 fZf p 



(22a) 



Here p = q— k, and c N (p) = \[2c p = v2c^. The above 
specified discretization of the wave numbers yields T q as 
polynomial 



F*[f x J] = [w x (q)/q 2 



kp 



fkfp 



(22b) 



One gets for the non-ergodicity parameters f q = <fi q (t — > 
oc) = (F q AA °°±F q AB °°)/w x (q) 

/ g x =^[A/]/{l + ^[A/]}- (23) 

There is no coupling between the fluctuations of the 
total density and the ones of the "charge" density. The 
mathematical structure of the two sets of equations for 
x = N and x = Z, respectively, is the same as the one 
studied previously for the density correlator <j> q (t) of a 
tagged particle in a simple liquid |B3fl, For the density 
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dynamics one also finds the small-g asymptote for the fre- 
quency Q^ 2 — {vTq) 2 + 0(q 4 ), reflecting free translation 
of the probe molecule. There also is the g _2 -divergency of 
the mode-coupling coefficients in , which implies the 
approach towards unity of the Lamb-Mossbauer factor 
for vanishing wave number: f^ = 1. For the "charge" 
dynamics, one gets a non-zero small-g limit for the fre- 
quency characterizing free rotation f^jj = 2v 2 { + 0(q 2 ). 
The mode-coupling coefficients do not diverge for q — > 0, 
since 6wz(q)/(Lq) 2 — > 1. Therefore, the non-ergodicity 
parameter for the variable pz(q,t) approaches a limit 



smaller than unity: f~_ 



< 1. 



C. The dipole correlator and the mean-squared 
displacements 

According to Eqs. (||), the knowledge of the dipole cor- 
relator C\ (t) and two of the mean-squared displacements 
Sr 2 (t) for a — A, B, or C is equivalent to the knowledge 
of the three independent elements of the symmetric ma- 
trix C(t). Using Eq. (0) and expanding Eq. (17a) for 
small wave numbers one gets 

d 2 C{t) - D + ftgC(i) + J / dt'm{t - t')d t ,C(t') = 0. 

Jo 

(24) 

This exact equation of motion for C(t) has to be 
solved with the initial condition from Eq. (^|). The fre- 
quency matrix is obtained as zero-wave-number limit 
from Eq. (p7b|): 



n 2 = (2v 2 R /L) 



Za 
zb 



-za 
-zb 



(25) 



Equation (f24|) implies C(0) 
one gets from Eq. (0): D = 



6 14 



2 uf, (z A + z B 



D + rjgC(O) = 0. Thus, 
6J + ftgCo, i.e., D 



ah 



The MCT approximation for the ke rnel m(i) is ob- 
tained by combining Eqs. (17t ) and (|l8b| ) and taking 
the zero- wave- vector limit. With Eq. flB17| ), one finds 



m(t) 



i 

6^ 



(26a) 

dfcfc 4 P S k clc h Jl h j k . (26b) 



Again, the theory simplifies considerably for symmetric 
molecules. In this case, one can transform Eq. ( J24) a s 
explained in connection with the derivation of Eq. ( 20J) . 
Using Eqs. ([l3|) one gets the exact equation of motion for 
the mean-squared displacement 



d?5rl(t) - 6«f + v\ / d£ m N (t - t') d v 8r 2 c {t') = 0, 
Jo 



(27) 



to be solved with the initial behavior Sr^it) = 3(i>t£) 2 + 
0(|£| 3 ). Similarly, one obtains for the dipole correlator 

d 2 C 1 (t) + 2v 2 R C 1 (t) 

+ 2vl I dt'm z (t-t')d t 'C 1 (t') = 0, (28) 
Jo 

to be solved with the initial decay C\{t) = 1 — (ur£) 2 + 
0(|i| 3 ). The MC T ap proximation for the kernels is ob- 
tained from Eq. (E6bh: 



(29a) 



i x (t) = F x [(f> x (t),<f>(t)}, x = NoxZ, 



T x \r\ f] = Qx / dk fc 4 pS k c JV (k) 2 w x (k) 31 fk, (29b) 
Jo 

where — l/{6n 2 ) and oez = L 2 /(72tt 2 ). Equa- 
tions (|2^) and (^) for the dipole correlator have the stan- 
dard form of the MCT equation. If <p^ (t) approaches zero 
for large times, the same approach towards equilibrium is 
exhibited by C\(t). If the solvent is a glass, f q > 0, and 
if the "charge" -density fluctuations (f>g (t) exhibit non- 
ergodic behavior, > 0, also the I — 1-reorientational 
correlator exhibits non-ergodic dynamics: 



d(i - co) = h = T z [f, f] I {1 + T z \f, /]}• (30) 

Parameter f\ is the long-wave-length limit of dis- 
cussed in Eq. @: f^ = ft. 

D. The quadrupole correlator 

The quadrupole correlator C 2 (i) = (3[e(t) • e] 2 - l)/2 
cannot be extracted from the correlators F^ b (t) with 
a, b = A or B. But let us consider a linear symmet- 
ric triatomic molecule. The third atom, labeled C, has 
its position in the center ?c- The preceding theory can 
be extended by adding as third variable the fluctuations 
p5 = exp(iq ■ fc)- The basic quantities are now the 
elements of the 3-by-3 matrix correlator, defined as in 
Eq. (||) with a, b — A, B or C. The correlator formed 
with pq{c[) = p~ + p~ — 2p S is a linear combination of 
the nine functions F^ h (t). An expansion for small q yields 



(PQm*PQ(9)) = m( < lL) 4 [C 2 (t) 



0(q% (31) 



In this case, C2{t) can be obtained in a similar manner 
as discussed above for C\(f). A diatomic molecule can 
be considered as a special mathematical limit of a tri- 
atomic one. Hence, there is no problem in principle to 
obtain C2(t) within a theory based on an interaction-site 
description. Motivated by this observation, an auxiliary 
site C shall be introduced |24|,|2j| and will be used as 
third basic variable. However, a complete theory with 
3-by-3 matrices shall not be developed. Rather some ad- 
ditional approximations will be introduced so that Cz(t) 
is obtained as corollary of the above formulated closed 
theory. 
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The quadrupole correlator can be written as small-g 
limit of a correlation functi on f ormed with tensor-density 
fluctuations defined in Eq. (Al) for q — (0, 0, q): C2(t) 



theory of Sec. [II C provides the results for the mean- 



lim g _ > o(P2(90) t)*P 2 (<lo))- Therefore, an exact Zwanzig- 
Mori equation can be derived as usual: 



d?C 2 (t)+6v 2 R C 2 (t) 



■6v£ 



dt'mf{t-t')d v C 2 {t') = 0. (32) 



The relaxation kernel m\^{t) is a correlator for fluctuating 
forces F^ 2 (q0, t) referring to angular-momentum index 
£ = 2 and helicity m = 0: 



lim(FR2(g0, t)*Jk2(g0)). 



(33) 



The time evolution of the fluctuating force is generated 
by the reduced Liouvillian C = QCQ, where Q projects 
perpendicular to p 2 (Qo) and Cp 2 (qo), and the Liouvillian 
C is defined by iCA(t) — dtA(t). The involved notation 
has been chosen in order to get the formulas in agreement 
with the ones of the more general theory in Ref. |[L3| . 
The procedure used for the theory of simple liquids |26j | 
shall be applied to derive an approximation for the kernel. 
First, the forces will be approximated by the projection 
onto the space of the simplest modes contributing F R2 — > 
V'F R2 . Here V projects onto the space spanned by the 
pair modes 



A a {k,p) = p^pp/y/NSp, a = A,BoTC. (34a) 

The essential step is the second one, where correlations 
of the pairs are replaced 

by products of correlations: (A a (k,p, t)*A a (k 1 ,p')) — > 
SUP$pP'{p%(t)*P%){pis(t)*pp)/NSp. This approxima- 
tion is done in particular for t = 0, thereby get- 
ting for the normalization matrix for the pair modes: 
(A a (k,p)*A b (k',p')} -> SjtpSpg'wf. Here w ab (k) = 
jo(k(z a — Zb)) generalizes Eq. (g) to a 3-by-3 matrix. As 
a result, one gets 



£p* abed 

The q — ► limit is carried out easily, reducing the sum 
over k and p to the one over k with k = —p. One obtains 
the kernel as mode-coupling functional 



\t) = 



dk Y, V£ b (k)Ff(t)<t> k (t)- (35) 

ab£{A,B,C} 



Let us restrict the discussion to symmetric molecules. 
For this cas e, an explicit expression for V 2 ab (k) is noted 
as Eqs. (A9) in Appendix |X| 

The correlators d>k(t) and F£ b (t) w ith a, b = A or B 
are taken from Sec. HI A| and Sec. Ill B| , respectively. The 



squared displacement Sr^if). The Gaussian approxima- 
tion shall be used to evaluate F^ c (t) w exp[— |fc 2 (5r^(t)]. 
The two remaining functions can be expressed in t erm s 
of tensor-density fluctuations according to Eq. (|A6|): 
^fe C W = EiVW+l)Mkza)MkO,t) for a = A or 
_B. As in the previous work E| , only the diagonal corre- 
lators shall be taken in the sum, i.e., the approximation 
(t)Kj (kz a )Fg G (t). 



will be used: F£ c 



IV. RESULTS 

A few concepts shall be mentioned which were intro- 
duced |2{J to describe the MCT-liquid-glass-transition 
dynamics. In the space of control parameters, a smooth 
function a is defined near the transition points, called 
the separation parameter. Glass states are characterized 
by a > 0, liquid states by a < 0, and a = defines the 
transition hypersurface. Suppose, only one control pa- 
rameter, say the density p, is varied near the transition 
point. Then one can write for small distance parame- 
ters e = (p — p c )/pc '■ c — Ce, C > 0. In addition to 
C, the transition point is characterized by a time scale 
to and by a number A, < A < 1. The scale to speci- 
fies properties of the transient dynamics, and A is called 
the exponent parameter. The latter determines a certain 
number B > 0, the critical exponent a, < a < 1/2, and 
the von-Schweidler exponent b, < b < 1. There are two 
critical time scales governing the bifurcation dynamics 
close to the transition: 



t /\a\ s , 4 = t B- 1/i /M 7 - 



(36) 



The anomalous exponents of the scales read: <5 = 
l/2o, 7 = l/2a+ 1/26. The hard-sphere-system (HSS) 
shall be used as solvent. There is only one control param- 
eter for the equilibrium structure, which shall be chosen 
as the packing fraction tp of the particles with diameter 
d, (p = irpd 3 /6. The distance parameter shall be given 
by the logarithm x of | e |: 



((p - cp c ) / tp c = ±10 



(37) 



The structure factor 5^ is calculated within the Percus- 



Yevick theory |20|. The wave numbers are discretized 
to M = 100 values with spacing hd = 0.4. For this 
solvent model, results for the density correlators and their 
spectra can be found in Ref. [E7||. The glassy dynamics 
is analyzed in Ref. |22|, from which one infers: tp c — 
0.516, C = 1.54, A = 0.735, a = 0.312, b = 0.583 and 
B = 0.836. Furthermore, t = 0.0236(d/i;) ||. 

Dumbbells of two fused hard spheres of diameters dA — 
ds = d shall be used as solute. The elongation parameter 
C = L/d quantifies the bond length. The solute-solvent- 
direct-correlation functions are also calculated within the 
Percus-Yevick theory. Within the tensor-density descrip- 
tion, the direct correlation functions ce(g) have been de- 
termined in Ref. |28|. These results are substituted in 
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the formulas of Appendix A, to evaluate the equilibrium 
structure in the site representation. In all summations 
over contributions due to various angular-momentum in- 
dices £, a cutoff £ co — 8 is chosen. It was checked for rep- 
resentative cases, that increasing the cutoff to £ co = 16 
does not significantly change the results to be discussed. 
The discretization of the various wave-vector integrals 
is done as sp ecifi ed ab ove for the solvent. The results 
in Sees. IV A and IV B deal wit h a sy mmetric dumbbell 
with via = tub = to, and in Sec. IV C , the molecule with 
ttia — 10m, mg = m is considered. 

Throughout the rest of this paper, the particle diam- 
eter is chosen as unit of length, d = 1, and the unit of 
time is chosen so that the thermal velocity of the solvent 
is v = 1. 



A. Structural arrest 

There are two control parameters for the system, 
namely the packing fraction ip of the solvent and the 
elongation £ of the solute molecule. Figure [l] exhibits 
the phase diagram. Phase I deals with states where 
ip is below the critical value ip c , i.e., the solvent is a 
liquid. In this case, the long-time limits of the mode- 
coupling kernels in Eqs. (|l8|) vanish. All solute correla- 
tors relax to zero for long times, and the molecule dif- 
fuses through the solvent. For ip > ip c , the solvent is 
a glass. Structural fluctuations behave nonergodically. 
In particular, a tagged solvent particle does not diffuse; 
rather it is localized. Since the atoms of the molecule 
with d,A = (1b = d experience the same interaction with 
the solvent as the solvent particles among each other, 
one expects the molecule to be localized as well. Indeed, 
Eq. ( p3| ) yields for ip > ip c : > 0. If ip increases from 



below ip c to above ip c , the long-time limit 



cNc 



oo) 

increases discontinuously at <p c from zero to > 0. 
Also, the quadrupole correlator exhibits nonergodic dy- 
namics: Ci(t — ► oo) = ji > 0. The cages surround- 
ing the molecule cause such strong steric hindrance, that 
quadrupole fluctuations of the orientational vector e can- 
not relax to zero. In this sense, the states ip > ip c are 
ideal glasses. 

There are two alternatives for the glass. Phase II deals 
with states for sufficiently small C- There is such small 
steric hindrance for a flip of the molecule's axis between 
the two energetically equivalent positions e and — e that 
Eq. (^) yields f^ = 0. The dynamics of the "charge" 
fluctuations is ergodic. In particular, the dipole correla- 
tor relaxes to zero: C\(t — > oo) = 0. Phase II is an amor- 
phous counterpart of a plastic crystal. For sufficiently 
large £, steric hindrance for dipole reorientations is so 
effective, that also the "charge" fluctuations behave non- 
ergodically. In this case, Eq. (23) yields a positive long- 
time limit < fg — cj)q (t — > oo). In particular, dipole- 
disturbances do not relax to zero: C\(t — > oo) = /i > 0. 
This phase III is a glass with all structural disturbances 



exhibiting nonergodic motion. The two phases II and III 
are separated by transitions at C, — Cc{^p), V > fc- With 
decreasing density, the steric hindrance for reorientations 
decreases. Thus, Cc(v) increases with decreasing (p, as 
shown by the full line in Fig. [l]. The transition curve 
terminates with horizontal slope at the largest critical 
elongation ( c = ((ip c ) = 0.380. Function ( c (ip) was calcu- 
lated before within the MCT based on the tensor-density 
description |2^| , and the transition curve of this theory is 
added as dashed line in Fig. [l| The results of the two the- 
ories are in qualitative agreement. It would be interesting 
if molecular-dynamics studies would decide, which of the 
two theories leads closer to reality. The asymptotic laws 
for the transition from phase II to phase III have earlier 
been described as type- A transition as can be inferred 
from Ref. ^] and the papers quoted there. At this tran- 
sition, C\(t — > oo) increases continuously with increasing 

The heavy full lines in Fig. ^| exhibit critical noner- 
godicity parameters c for £ — 0.80, calculated from 
Eq. (E3I) for the liquid-glass transition point ip = tp c . 
These quantities are Lamb-Mossbauer factors of the 
molecule. The function f^ c can be measured, in prin- 
ciple, as cross section for incoherent neutron scattering 
from the solute, provided both centers A and B are 
identical atoms without spin. As expected for a local- 
ized probability-distribution Fourier transform, the c - 
versus-g curves decrease with increasing q. Most remark- 
able are the kinks exhibited by f^ c for wave numbers q 
near 5, 12.5 and 20, and by f^ c for q near 10 and 17.5. 
The light full lines exhibit c calculated with Eq. ( |A7| ) 
from the critical nonergodicity parameters f c (q£0) J13| . 
The results of both approximation theories are in semi- 
quantitative agreement, in particular concerning the po- 
sition and size of the kinks. The / c (gl0)-versus-g curves 
are bell shaped, close to Gaussians p3[ . They enter 
Eq. (A?) with prefactors jg(q(/2) 2 = 0{q 2e ), so that the 
maximum of the contribution from angular-momentum- 
index £ occurs at some g™ ax which increases with £. The 
separate contributions for different £ are shown as dotted 
lines in Fig. ^. Thus, the kinks are due to interference 
effects of the f c (q£0) with the intramolecular form fac- 
tors ji(q(/2). Let us add, that also the Lamb-Mossbauer 
factors of the atoms, /° c , are well described by Gaus- 
sians for q < 10; in particular these functions do not 
exhibit kinks. Figure || demonstrates for a case of strong 
steric hindrance for reoricntational motion that angular- 
momentum variables for £ up to 6 are relevant to describe 
the arrested structure, and that the description of the 
molecule by site-density fluctuations properly accounts 
for the contributions with £ > 2. 

Figure || exhibits c representative for weak steric 
hindrance for the reorientational dynamics. Naturally, 
the contributions due to the arrest of fluctuations of ten- 
sor densities with large £ are suppressed. The contribu- 
tions for £ — and £ = 2 are sufficient to explain f^ c , in 
particular its kink for q near 12.5. Similarly, the contri- 
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tuitions for I = 1 and £ = 3 are necessary and sufficient to 
explain f^ c with its kink for q near 17.5. The dynamics is 
strongly influenced by precursor phenomena of the tran- 
sition from phase II to phase III. This is demonstrated, 
for example, by the strong decrease of f? = fg^ for the 
result shown in the lower panel of Fig. |3j in comparison 
to the one shown in the lower panel of Fig. ||. The two 
approximation theories under discussion yield different 
numbers for the value £ c for the transition point. It is 
meaningless to compare different approximations for re- 
sults near a singularity Cc, referring to the same value 
It is more meaningful to compare results for the same 
relative distance from the critical point, — £c)/Ccj as 
is done in Fig. [| Let us mention that f^ c shown by the 
heavy and light full lines would be somewhat closer, if 
one had compared elongations yielding the same value 
for/f. 

Figure || exhibits critical Lamb-Mossbauer factors f£ c 
as function of the elongation. The lower panel demon- 
strates the transition from phase II for £ < ( c to phase 
III for £ > ( c . For strong steric hindrance, say £ > 0.8, 



fq C is rather close to f Zc provided q is not too small, 
say q > 3. For ( approaching £ c , the f^ c fall below 
fq C . Most remarkable are the wiggles or even minima 
of the curves. These are the analogues of the kinks, dis- 
cussed above in connection with Figs. || and ||. Again, 
these anomalies can be explained as interference effects 
between the geometric structure factors je(q(/2) 2 and the 
nonergodicity parameters f c (q£0) according to Eq. (A7). 



Let us consider for an intermediate wave vector as 
shown for curves b or c in the upper panel of Fig. |J. For 
small C, say C_< 0.4, the £ = contribution dominates the 
sum in Eq. (A7), as can be inferred from Fig. ||. Func- 
tion f c (q00) reflects the isotropic part of the arrested 
fluctuations, and hence it is practically ^-independent 
(as shown in Fig. 5 of Ref. (If]). Since jo(gC/2) 2 = 
1 ~ TlilO 2 + 0((q() 4 ) decreases with increasing the 
/ ( j Vc_versus_ C curve decreases too; and the decrease is 
stronger for larger q. The / c (g20) increase from for 
( = to values near 0.5 for £ = 1 (as shown in Fig. 5 of 
Ref. jl3)). Also, the geometric structure factor increases 
strongly with C : j 2 (<zC/2) 2 = ((gC) 2 /60) 2 + O((<z0 6 )- 
The combined effect of both increases causes the increase 
of the / ( j Vc -versus-C curve for larger (. The resulting 
minimum occurs for smaller £ if q is larger, and this ex- 
plains the difference between the two curves b and c. The 
theory produces the minima since it accounts for the ar- 
rest of tensor- density fluctuations for £ > 2. 



B. Correlation functions and spectra near the glass 
transition 

Figure |5| demonstrates the evolution of the dynam- 
ics for the correlators <j>^ (t) and <j>^ (t) for intermediate 
wave numbers q near the transition from phase I to phase 
III. The oscillatory transient dynamics occurs within the 



short-time window t < 1. The control-parameter sensi- 
tive glassy dynamics occurs for longer times for packing 
fractions ip near if c . At the transition point <f = tf c , 
the correlators decrease in a stretched manner towards 
the plateau values f*° as shown by the dotted lines. In- 
creasing if above ip c , the long-time limits increase, as 
shown for (f>g (t) for q = 7.4 and C = 0.80. Decreasing ip 
below if c , the correlators cross the plateau at some time 
Tj3, and then decay towards zero. The decay from the 
plateau c to zero is the a-process for (t) . It is char- 
acterized by a time scale r a , which can be defined, e.g., 
by 4>q{T a ) — fg C /2. Upon decreasing <f c — ip towards 
zero, the time scales Tp and r a increase towards infinity 
proportional to t a and t' a , respectively, cited in Eq. (|3^). 
The figure exemplifies the standard MCT-bifurcation sce- 
nario. For small | ip — f c |, the results can be described in 
terms of scaling laws. This was explained in Refs. |p2 23 
for the HSS, and the discussion shall not be repeated 
here. 

One can deduce from Fig. || that for £ = 0.80 and 
q > 5 the plateaus for both types of density fluctuations 
are very close to each other: f^ c w f^ c . The upper two 
panels of Fig. || demonstrate that also the dynamics is 
nearly the same, <t>q{t) ~ 4>q(t)- This means that for 
q( > 4 and for strong steric hindrance, the cross cor- 
relations F^ B (t) are very small. The reason is that the 
intramolecular correlation factors ji{qC/^) are small, and 
thus interference effects between the density fluctuations 
of the two interaction sites are suppressed. Coherence ef- 
fects can be expected only for smaller wave vectors. For 
this case, the functions can be understood in terms of 
their small-q asymptotes, Eq. (0). 

The lower two panels in Fig. ^ deal with weak steric 
hindrance. In this case, the "charge" -density fluctuations 
behave quite differently from the number-density fluctua- 
tions. The most important origin of this difference is the 
reduction of the mode-coupling vertices V q z kp relative to 



V^ kp in Eq. ( |22b| ) . For small elongations of the molecule, 



the effective solute-solvent potentials for reorientations 
are small. Therefore, the f Zc decrease strongly relative 
to fq C for C decreasing towards Co a s is shown in Fig. ||. 
Upon approaching ( c , the a-peak strength of <j>g(t) gets 
suppressed relative to the one of <f) q (t) . Within phase II, 
the "charge" -density fluctuations relax to zero as in a nor- 
mal liquid. This implies as precursor phenomenon, that 
the time scale t z of the "charge" -density-fluctuation a- 
process shortens relative to the scale for the number- 
density fluctuations. Thus, the small-£ behavior shown 
in the lower two panels of Fig. || is due to disturbances 
of the standard MCT-transition scenario by the nearby 
type-A transition. 

The correlators C\(t) and C^if) are shown in Fig. ^ 
for the critical point ip = f c and for two liquid states 
near the transition from phase I to phase III. For £ = 
0.80, the anisotropic distribution of the solvent particles 
around the molecule leads to a stronger coupling to the 
dipole reorientations than to the reorientations for the 



8 



quadrupole, and therefore the plateau for the former is 
higher than for the latter, /f > f%. A leading order ex- 
pansion of the solutions of the equations of motion ( |2§| ) 
and ( [32|) in terms of the small parameter Ce(t) — ff leads 
to the factorization in the critical amplitude hf, and a 
function G(t) called the /3-correlator, C^{t) — f? — hgG(t). 
The latter is the same for all correlation functions. It 
obeys the first scaling law of MCT, G(t) = ■ s /\a\ g±(t/t a ) 
for cr^O. The master functions g±(t) are determined by 
A. They also describe the dynamics of the solvent in the 
window where | 4> q (t) — f° |<S 1 |2l| . In particular there 
holds g_(i_) = 0, i_ = 0.704, so that both correlators 
Ci(t) cross their plateau at the same time Tp = t^t a . The 
non-linear mode-coupling effects request, that the corre- 
lators approach zero roughly at the same time. Thus 
one understands the general differences between the a- 
processes, which were mentioned in the introduction: 
the a-process for dipole relaxation is stronger, slower, 
and less stretched than the one for quadrupole relax- 
ation. This finding is in qualitative agreement with the 
ones of the theory based on tensor-density representa- 
tion of the structure fl4| . There are, however, quantita- 
tive differences between the two approximation schemes. 
The plateaus /f = 0.905 and f% = 0.674 are smaller 
than the corresponding values 0.943 and 0.835 found in 
Ref. Q and the amplitudes hi = 0.19 and h 2 = 0.40 
are bigger than the corresponding values 0.13 and 0.35 
calculated previously juj. The times T l a characteriz- 
ing the a-process shall be defined by Gi>(t^) — ff/2. 
They are marked by open squares in Fig. 0. The val- 
ues tI = 5.21 x 10 5 , tI = 1.64 x 10 5 for x = 3 and 
C = 0.80 are smaller than the ones reported in Ref. p| . 
The present theory implies a somewhat weaker coupling 
of the reorientational degrees of freedom of the molecule 
to the dynamics of the solvent, than found earlier |jl4) . 
This holds also for the small elongation £ = 0.43. The 
approach towards the transition from phase III to phase 
II, leads to a suppression of /f , as discussed for the f^ c 
in Fig. ||. The dipole relaxation speeds up for £ — * £ c , as 
discussed for the lower panels of Fig. [|. This is reflected 
by an enhancement of h\ — 1.60 relative to the ampli- 
tudes cited for ( = 0.80 but also relative to h 2 = 0.49. 

One can perform lim^o- i'qitt'o) — <f>q(t) for the so- 
lutions of Eq. (|l4|), where <j> q {t) can be evaluated from 
the mode-coupling functional at the critical point. It 
obeys the initial condition (f> q (t) = /| — h q t b + 0(t 2b ). 

Function (j) q (t) can be considered as shape function of 
the a-process, and the result implies the second scaling 
law of MCT, also referred to as superposition principle: 
4> q (t) = <t> q (t/t' a ) for — > — ^(|. Corresponding laws 
hold for all functions, as is demonstrated in detail for 
the HSS in Refs. p2] , p3| . In particular, one gets for the 
reorientational correlators for a — > 0— : 

C e (t)=C e (t), t = t/t' (T , t a <t, (38a) 

and this corresponds to the a-scaling law for the suscep- 
tibility spectra 



Xi'(w) = xi'(w), w = u<., w<l/V (38b) 

The initial decay of the master function Cg(t) is described 
by von Schweidler's law 

Ci{t) = f c e - h t t\ t^0, (39a) 

which is equivalent to a power-law tail of the master spec- 
trum x"{£u): 

X'l(u>) = h e sin(7r&/2) T(l + b)/Cj b , Co 00. (39b) 

The upper panel of Fig. ^ exhibits the a-process mas- 
ter spectra for the reorientational processes for Q = 0.80 
and for the dimensionless longitudinal elastic modulus 
m q= o(t) of the solvent. The latter can be measured 
by Brillouin-scattering spectroscopy. It probes a tensor- 
density fluctuation for i = 0. The von-Schwcidler-law 
tails describe the spectra for frequencies exceeding the 
position a) max of the susceptibility maximum by a factor 
of about 100, as shown by the dashed lines. Since > /| 
and both plateau values are rather large, one underst and s 
from the theory for the leading corrections to Eq. ( |39b| ) 
p2[ that for decreasing Co the von-Schweidler asymptote 
underestimates the spectrum, and this by larger values 
for £ = 1 than for I = 2. For smaller plateau values, the 
von-Schweidler asymptote may overestimate the spec- 
trum, as is exemplified for the modulus. The lower panel 
of Fig. ^ demonstrates that the a-processes speed up if 
steric hindrance is decreased. As precursor of the tran- 
sition to phase II, the spectrum for the dipole response 
is located at much higher frequencies than that for the 
quadrupole response. Traditionally, dielectric-loss spec- 
tra have been fitted by the ones of the Kohlrausch-law 
<f>K(t) — ^4exp[— {tB)P\. Such fits also describe a ma- 
jor part of the spectra in Fig. |?], as shown by the dotted 
lines. The parameters A and B are adjusted to match 
the susceptibility maximum. The stretching exponent (3 
is chosen so that the spectrum is fitted at half maximum 
x"((2) max )/2. If one denotes the width in log 10 w at half 
maximum by W, stretching means that this parameter is 
larger than the value Wn =1.14 characterizing a Debye 
process, </>z>(i) = exp(— t). The upper panel of Fig. 
quantifies the general results of the theory for strong 
steric hindrance: Xi(^max) > ^'("D. 
and Pi > [32- It quantifies also the forth property cited 
in the introduction <I^ lax < liJJ,^. A further general prop- 
erty is @ 2 > Po- 

The discussion of power-law spectra is done more con- 
veniently in a double logarithmic diagram as shown in 
Fig. |^ for normalized dipole-fluctuation-a-process spec- 
tra C{'(d))d> max /7 1 c = xi(^)^max//i c w as function of 
Co/Cdmax- One notices, that there is a white-noise spec- 
trum for Co below <D max . The high-frequency wing of the 
Kohlrausch-law fit decreases proportional to Co~@ and 
underestimates the spectrum Xi(^) considerably. Be- 
cause of the von-Schweidler asymptote, which is shown as 
dashed straight line, the spectrum exhibits an enhanced 







high-frequency wing. Dixon et al. |3l| made the remark- 
able observation that their dielectric spectra could be col- 
lapsed on one master curve, if the vertical axis is rescaled 
by w^ 1 and the horizontal axis by + w _1 ), where 

w = W/Wd- In Fig. N this scaling is used and the data 
for glycerol from Ref. plj] are included. The spectra for 
molecules with £ = 0.80 and ( = 0.60, which are relevant 
for the description of van der Waals systems [jl4| , follow 
the mentioned scaling surprisingly well. This finding ap- 
pears non-trivial, since the scaling is not reproduced by 
the MCT results of the basic quantities <j> q (t) @. The 
rescaled spectrum for £ = 0.43 deviates from the scaling 
law for + w^ 1 ) log 10 (u)/u) max ) > 5. 

It might appear problematic that the dipole correlator 
was calculated within a different approximation scheme 
than the quadrupole correlator. But, there is no difficulty 
to a lso ev aluate C\{t) within the scheme explained in 
Sec. [II D for the evaluation of C2(t). Figure ^ presents 
a comparison of C\(t) obtained along the two specified 
routes. The two results for the small elongation ( = 0.43 
are close to each other. The discrepancies are mainly 
due to the 7% difference between the two plateau values 
fi- With increasing £, the discrepancies decrease. For 
the large elongation £ = 0.80, the results are practically 
undistinguishable . 



C. Structural relaxation versus transient dynamics 

Let us introduce Fourier-Laplace transforms of func- 
tions of time, say f(t), to functions of frequency, say 
f(u>), with the convention f(ui) = i dtexp(izt) f(t), 
z = u> + iO. The equations of motion ( |l7j ) with the initial 
conditions from Eq. (H) arc transformed to 



Similar reasoning leads from Eq. (HJ) to 



fl 2 „ ra q {uj)\ [u F q (iu) + w 9 ] - fl 2 „ F,(w) = 0. (40) 



Within the glass, the long-time limits of F q (t) and m q (t) 
do not vanish, i.e., the transformed quantities exhibit 
zero- frequency poles. One gets, for example, u>F q (u>) — > 
— for uj — ► 0, where the strength — F^° of the poles 
follow from Eq. ([l9]). Continuity of the solutions of the 
MCT-equation of motion implies that for vanishing fre- 
quencies and for vanishing distances from the transition 
points, xn. q (u>) becomes arbitrarily large. Hence, combi- 
nations like ui + i^ q with constants £ q can be neglected 
compared to mJu). Therefore, in the region of glassy 
dynamics, Eq. ( |40|) can be modified to 



F q (u)) - m q (u)vf q =i£ +umq(u) F q (uj). 



(41) 



Let us assume, that this equation has a solution, to be 
denoted by F*(w), which is defined for all frequencies, 
so that it can be back transformed to a function F*(t), 
defined for all £ > 0. Choosing £ q properly, Eq. ( pLt| ) can 
be written as 



F; (t) - m* (t) w, = - (d/dt) [ dt' m* (t - f ) F* (f 

Jo 



P q (t)-m q (t) = -(d/dt) dt'm* q {t-t')p q {t'). (42b) 
Jo 

These formulas have to be supplemented with the MCT 
expressions for the kernels 

m* (t) = T q [F* (i) , 0* (t)] , m* (i) - T q W (t)\ . (42c) 



Equations (|4j) for the glassy dynamics are scale invari- 
ant. With one set of solutions 4>* q (t) and F*(i), also the 
set <j) q u (t) = <j>* q (ut) and F*"(i) = F*(ut) provides a solu- 
tion for arbitrary u > 0. To fix the solution uniquely, one 
can introduce positive numbers y q and positive definite 
matrices y q to specify the initial condition as power-law 
asymptotes Jl(| 

p;W (^Ao) 1/3 - y q , <t>* q (t) (^Ao) 1/3 - y q , (t/*o) - o. 

(43) 

The theory of the asymptotic solution of the MCT 
equations for simple systems had been built on the ana- 
logue of Eq. ( pi] ) with £g neglected [p2| . The present 
theory extends the previous one by the introduction of 
matrices. It seems obvious that the previous results for 
asymptotic expansions hold in a properly extended ver- 
sion. Let us only note the formula for the long-time decay 
of the correlators at the critical point <p — ip c p3,E3|: 



Fg(t) 



TTlOO C 



• U q (t /t) a [l + Kq(t Q /t) a + O((t /t) 2a )}. 

(44) 



The exponent a, mentioned above, can be calculated from 
the mode-coupling functional at the critical point. The 
same is true for the plateau values F^° c , and the ampli- 
tudes H g and K q . The dependence of the solution from 
the transient dynamics is given by the single number to- 
Let us anticipate that Eq. ( |44| ) and similar results can be 
extended to a complete solution. One concludes that the 
glassy dynamics is determined, u p to a scale to, by the 
mode-coupling func tionals in E q. ( 42c ). 

Equations (Bll) and (B17) show, that the mode- 



(42a) 



coupling functionals T q and T q b are specified by the den- 
sity p, the static structure factor S q , the direct correla- 
tion function c q of the solvent, and the solute-solvent 
direct correlation functions c ql i.e., by equilibrium quan- 
tities. They are the same for systems with a Newto- 
nian dynamics, as considered in this paper, and for a 
model with a Brownian dynamics, as is to be used for 
the description of colloidal suspensions. In particular, the 
mode-coupling functionals are independent of the parti- 
cle masses to, tua and ms- Thus, the glassy dynamics of 
the molecule in the simple liquid does not depend on the 
inertia parameters specifying the microscopic equations 
of motion. The same conclusions on the glassy dynam- 
ics, which were cited in the introduction for the basic 
version of MCT, hold for the model studied in this pa- 
per. Let us add, that neither the temperature T, nor 
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the interparticle-interaction potentials V of the solvent, 
nor the solute-solvent-interaction potentials V a enter ex- 
plicitly the mode-coupling functionals. These quantities 
only enter implicitly via S q , c q , and . 

The independence of the glassy dynamics on the iner- 
tia parameters is demonstrated in Fig. |l^ for four states 
of the liquid. It is shown that the reorientational dy- 
namics of the dipole does not change for t > 1 even if 
the mass ratio of the atoms m a / 'mg is altered by a fac- 
tor 10. The transient dynamics, which deals with over- 
damped librations, exhibits an isotope effect. There is 
no fitting parameter involved in the shown diagram. The 
scale to neither depends on the density of the solvent nor 
on the elongation parameter, and Eqs. (|4|) describe the 
complete control-parameter dependence of the glassy dy- 
namics. 



V. SUMMARY 

The MCT for simple systems with a dilute solute of 
atoms has been generalized to the one with a dilute so- 
lute of diatomic molecules. The derived equations of mo- 
tion generalize the ones for atoms in the sense that scalar 
functions are replaced by 2-by-2 matrix functions. These 
generalizations result from the description of the position 
of the molecule in terms of interaction-site-density fluc- 
tuations. The numerical effort required for a solution of 
the equations is not seriously larger than the one needed 
to solve the corresponding equations for atomic solutes. 
This holds in particular for symmetric molecules where 
the matrix equations can be diagonalized by a linear 
transformation to number-density and "charge" -density 
fluctuations. The theory implies that the dynamics out- 
side the transient regime is determined, up to an overall 
time scale to, by the equilibrium structure. In particular, 
it is independent of the inertia parameters of the system. 

It was shown that the present theory reproduces all 
qualitative results obtained within the preceding much 
more involved theory based on a description of the so- 
lute by tensor-density fluctuations. Both theories yield 
a similar phase diagram, Fig. ^. The characteristic dif- 
ferences of the reorientational correlators exhibited for 
strong steric hindrance of rotations as opposed to weak 
steric hindrance are obtained here, Fig. g, as previously. 
There are systematic quantitative differences between the 
two approximation approaches in the sense, that the im- 
plications of the cage effect are weaker in the present 
theory than in the preceding work The earlier 

findings on the differences of the spectra referring to re- 
sponses with angular-momentum index £ = 0, 1 and 2 
have been corroborated by separating the a-peaks from 
the complete spectra with the aid of the a-scaling law, 
Fig.@. 

The critical non-ergodicity parameters f^ c have been 
calculated, which determine the form factors of quasi- 
elastic scattering from the liquid near the glass transition. 



Contrary to what one finds for atomic solutes, these are 
not bell-shaped functions of wave numbers q, rather they 
exhibit kinks, Fig. 0. The form factors for some wave 
numbers q vary non-monotonically with changes of the 
elongation of the molecule, Fig. [|. Analyzing these find- 
ings in terms of form factors defined for fixed angular- 
momentum index £, one can explain the results as due 
to intra-molecular interference effects demonstrating that 
the theory accounts for reorientational correlations with 
angular-momentum index £ > 2. 

It was shown that the a-relaxation spectra for dipole 
reorientations with strong steric hindrance obey the scal- 
ing law proposed by Dixon et al. pl[ | within the window 
and within the accuracy level considered by these au- 
thors. There is no fitting parameter involved in the con- 
struction of Fig. ||, which demonstrates this finding for 
the two elongation parameters £ = 0.60 and 0.80. Thus, 
it is not justified to use the cited empirical scaling as an 
argument against the applicability of MCT for a discus- 
sion of dielectric-loss spectra. 
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APPENDIX A: TENSOR-DENSITY 
REPRESENTATIONS 

Following the conventions of Refs. |lj| and |L4|], nor- 
malized tensor-density fluctuations shall be used by de- 
composing the molecule's position variable in plane waves 
exp(iq-fc) for the center of mass fc and in spherical har- 
monics Y f m (e ) for the orientation vector e: 

pT (q) = i l yfa exp(iq ■ f c ) Y e m (e). ( Al ) 

The molecule-solvent interactions are specified by a set 
of direct correlation functions 

Ci(q) = {p% P°M)) I (pS q ); % = (0, 0, g). (A2) 

The structure dynamics is described by the matrix of 
correlators, defined by 

<f>tk(qm,t) = (pT(<b,t)*p%(qo)). (A3) 

Since the position vectors of the interaction sites can 
be written as f a — fc + z a e, the Rayleigh expansion of 
the exponential in Eq. (|l|) yields the formula 
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(A4) 



Substitution of this expression into Eq. (^|) leads with 
Eq. (A2) to an expression connecting the direct correla- 
tion functions c a q with Ci(q): 



E < i = E V2l+Tj e (qz a ) ce(q). 



(A5) 



Substitution of Eq. ([A|) into Eq. (|) leads with Eq. ([Af) 
to the expression of the site-density correlators in terms 
of the tensor-density correlators: 

= E VW+TWkTT) Ji(qz a ) j k (qz b ) fa k (qO,t). 



Ck 



(A6) 



For the long-time limits of the correlators Fg (t) one 
gets a combination of the non-crgodicity parameters 
fikiqm) — (f>ek(qfn,t — > oo). If one uses the diagonal 
approximation fik(q,m) = dikf(qim), one can write 

Ff(t -> oo) = Y^W + l)h(qZa)Mqz b ) f(q£0). (A7) 

t 

Su bstit uting the Rayleigh expansion of Eq. (Q) into 
Eq. (34a), one can express the pair modes in terms of 
those formed with tensor-density fluctuations: 

A a (k,p) = ^E^^) 1 ™ [pT(k) P f/VNS~ P ). 



I'm 



(A8) 



Therefore, the overlaps of the forces with the pair modes 
can be expressed as sums of the corresponding quantities 
calculated in Ref. [Ol. One finds for the mode-coupling 
coefficients in Eq. (|35|): 



V 2 ab (k) = [fc 2 P 5 fc /2^]E(w,7 1 ) ac ^c(fc) J D d (fc) (w-y b 

cd 

(A9a) 

Here one gets for a = A, B or C 

D *( k ) = T2 E(- 1 )" (£+J) ( 2£ + !) V2J + lj e (kz a ) cj(k) 



u 



x [J(J + I) + 6 -£(£+1)' 



2 £ J 




, (A9b) 



in Ref. |18|] by extending the procedure used originally 
for atomic systems ^(|. But the reported formulas JlJ] 
do not seem appropriate, since they do not reduce to the 
ones for tagged particle motion if the limit of a vanishing 
elongation parameter £ is considered. Therefore, an al- 
ternative derivation will be presented which starts from 
the theory developed by Mori and Fujisaka |3j| for an 
approximate treatment of nonlinear fluctuations. The 
application of this theory will be explained by rederiv- 
ing the equations for the solvent before the ones for the 
solute are worked out. 



1. The Mori-Fujisaka equations 

Let us consider a set of distinguished dynamical vari- 
ables A a ,a = 1, 2, defined as functions on the sys- 
tem's phase space. The time evolution is generated by 
the Liouvillian, A a (t) = exp(iCt)A a . Using the canon- 
ical averaging to define scalar products in the space of 
variables, (A, B) = (A*B), the Liouvillian is Hermi- 
tian. It is the goal to derive equations of motion for 
the matrix of correlators C a p(t) = (A a (t) : Ap). The 
initial condition is given by the positive definite ma- 
trix g a p = (A a ,Ap). The theory starts with a gener- 
alized Fokker-Planck equation for the distribution func- 
tion g a = H a S(A a — a a ), a — {a\, a 2 , ■ ■ •)■ It is assumed 
that the time scales for the fluctuations of the A a and 
their products is larger than the ones for the Langevin 
fluctuating forces. The spectra of the latter can then be 
approximated by a constant matrix T a p. It is assumed 
furthermore, that the T a p are independent of the distin- 
guished variables, and that the equilibrium distribution 
of the latter is Gaussian: 

(9a) = C expf-^E a *9ap a *p)- ( B1 ) 

In the cited Fokker-Planck equation there occurs the 
streaming velocity v a (a) given by 

v a (a){g a ) = (A a9a ) =k B Tj2-^-({^M9a). (B2) 

The Fokker-Planck equation is now reduced by project- 
ing out the subspace of distinguished variables. There 
appears the frequency matrix specifying the linear con- 
tribution to the streaming term, 



where the last factor denotes Wigner's 3-j symbol. 



n al3 =^2(A a ,jCA y )g-^. 



(B3) 



APPENDIX B: MODE-COUPLING 
COEFFICIENTS 

Mode-coupling equations based on a description of the 
molecules by site-density fluctuations have been derived 



The nonlinear contributions enter as combination 
fa = J dav a (a) g a - i E ^apAp . 



(B4) 







They determine the relaxation kernel as 
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M a/3 (t) = £(/«(*), A) s7/s- 



(B5) 



The time evolution is generated by the reduced Liouvil- 
lian CJ — QCQ, f a (t) — exp(iC't)f a , where Q is the pro- 
jector on the space perpendicular to the one spanned by 
the distinguished variables. The result, which is equiva- 
lent to Eq. (2.15) in Ref. f§, reads 



dtC afi (t) = -^[(iQ a7 + r Q7 )c 7/3 (t) 



+ / dt' M Q1 {t-t')C-yp(t')}. (B6) 
Jo 



2. MCT equations for simple systems 

To get a description of slowly varying structural 
fluctuations, the original reasoning of MCT shall be 
adopted, and as distinguished variables the density fluc- 
tuations and the longitudinal current fluctuations 
jq- = J2 K v k,z exp(iq ■ r K ) will be chosen. The variable 
label of the preceding subsection consists of two bits, 
a = (X,q), A = 1,2, so that A\$ = p$, A 2q - = jq. Nota- 
tions follow the ones of the first paragraph of Sec. 0. 
One gets g lq - lq ~ = NS q , g 2 q2q = Nv 2 , VL lq - 2 q = q, 
^2qiq = ^Iq/q, an d all the other elements of the matrices 
g and fl are zero. Because of translational invariance, 
one gets C^. = unless q — k, and rotational in- 
variance implies that C\^^{t) depends on the modulus q 
only. Th e same holds for the matrices T and M, so that 
Eq. (B6) reduces to a two-by-two matrix equation with 
q appearing as a parameter. Since Cpq* is an element of 
the distinguished set of variables, the kernels T a p and 
M a p(t) vanish unless a = f3 = (2, q). The latter shall be 
denoted by T q and M q (t), respectively. Equation (jB6| ) 
can then be reduced to the equation of motion for the 
normalized density correlator <fi q (t) = Ciq-iq-(t)/NS q : 

d^ q (t) + T q d t 4> q (t) 

+ fl 2 q <j)q(t) + [ dt' Mq(t - t') d t '<p q (t') = 0. (B7) 

Jo 

The application of the Mori-Fujisaka formalism can be 
summarized as follows. The relaxation kernel fl 2 m q (t) 
of the exact Eq. ( |l4| ) is approximately split into a white 
noise contribution 2T q 8(t) and a remainder M q (t), where 
well defined formulas for the two contributions are avail- 
able. In this paper, the kernel 



Eq. (B5) for the kernel, 



T q is neglected and 



Mq(t) = (ftft) , fq) /NV 2 , ff = fa, 



(B8) 



shall be approximated further. 

If one writes = pq and a 2q - = jq-, one can denote 
the equilibrium distribution w{a) = (g a ) for the solvent 
variables as 



;(o) = C exp|-(l/2iV) £[(1 - pc q ) p 9 p 9 * 



+ (l/v 2 )kfq]\- (B9) 



This is used to derive from Eqs. (B2) and (B4): 

f 9 = -i(pv 2 /N) ]T(fc ■ q/q) c k p % p ? _ % + Sfq-, (BIO) 



where Sfg denotes a term whose contribution to the 
memory kernel turns out to be irrelevant for the struc- 
tural relaxation processes, and therefore shall be ne- 
glected. The remaining task is the evaluation of averages 
{Pk{t)pp(t)> Pk'Pp')' wnere the time evolution is generated 
by the reduced Liouvillian. Here the original MCT ansatz 
is used: (pt(t)pp){pp(t)ppi) + (k «-> p). As a result, one 
gets M q (t) — £l 2 m q (t) with the well known expression 
for the mode-coupling functional in Eqs. (II 



T q [f] = (/V16ttV) / dkSqS k S p 



[q- kc k + q- pep} 2 fa f p , (Bll) 



with p = q — k. 



3. MCT equations for the solute molecule 

It is straightforward to generalize the preceding deriva- 
tion to systems with a diatomic solute molecule provided 
the latter is considered as flexible. Therefore, let us use 
this modification of the problem. It will be assumed that 
the kinetic energy of the molecule is Y> a (rn a /2)v 2 and 
that there is a binding potential V^{\ rA — fs between 
the two interaction sites. The exact Eq. (17a) remains 



flexible molecule 2 
the simple velocity corr elat or J' q ' 



valid with the frequency matrix replaced by t hat of the 
It is defined via Eq. jl7b| ) with 
b = 6 ab (k B T/m a ). 
The formulas of Sec. [B 1| shall be applied with the ex- 
tended set of densities (p<f, pf,p^) and the corresponding 
longitudinal currents (jf, j~ , j§). Thus the index con- 
sists of three bits a = (r,X,q), where A = 1,2 discrim- 
inates between densities and currents, and r = 0,A,B 
indicates solvent, atom A and atom B, respectively. In 
the infinite dilution limit N — > oo, the equations for the 
solute do not directly couple to those for the solvent. 
Therefore, Eq. (B7) remains valid and one gets a modi- 
fication of Eq. (17a) for the solute: 



d 2 Fq(t)+rqd t Fq(t)+nf q e2 F q (t) 

+ f dt'M q (t-t')d t >F q (t') = 0. (B12) 
Jo 



The relaxation kernel reads 



Mf(t) = (/£(*), fi) (m b /k B T), /f = f a2q -. (B13) 
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The determination of the Gaussian distribution of the 
extended distinguished variables requires the inversion 
of the three-by-three matrix (A T \f, A a iq)- Making use of 
the infinite dilution limit N — > oo, one gets 

(g a )=w°(a) exp|(l/2)^^[(2p/Ar)^ c «p f pf 

q ab 

- (w- 1 )^ - (m a /k B T) 6 ab ~j$ jf ]}, (B14) 

where w (a) is the distribution for the solvent variables 
given by Eq. (BE). This expression is used to work 
out t he fl uctuating force as explained in connection with 
Eq. O: 



fS = -i{p/N){k B T/m a ) y\(q - k) ■ q/q] c a k pip.. 



k 



k ^q-k 



(B15) 



As above, the contributions due to <5/2 are neglected and 
the remaining pair correlations are factorized. This leads 
to 



M q (t)=n( e2 vn q (t), 



(B16) 



where the functional for the kernel in Eq. (18a) reads 



ffll f] = <r 2 E < (p/ 8 ^ / * 



x s. 



Jkfp, (B17) 



with p abb revia ting q — k. Neglecting the friction term 
T q in Eq. the MCT equations (||) and © are 

derived for the flexible molecule. 

Obviously, a theory for a rigid molecule can be ob- 
tained from one for a flexible molecule only within a 
quantum-mechanical approach. One has to consider the 
case where excitation energies for translational and rota- 
tional motion are small compared to the thermal energy, 
while the energies for vibrational excitations are large. 
Let us assume, that the formulas can be obtained by 
replacing all equilibrium averages in the preceding equa- 
tions by the correct quantum mechanical ones, where the 
latter can be evaluated for the classical molecule model 
with five degrees of freedom. This amounts to replacing 
structure functions by the classical quantities, in partic- 
ular the replacement of Ct^ e by H q . Thereby Eq. (B12) 
produces Eqs. (|i~7|). 
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FIG. 1. Phase diagram of a dilute solute of symmetric 
dumbbell molecules with elongation £ consisting of two fused 
hard spheres which are immersed in a hard-sphere system 
(HSS) with packing fraction <p. The horizontal line marks 
the liquid-glass transition at the critical packing fraction 
fc = 0.516. The other full line is the curve C,c(tp) of criti- 
cal elongations for a type-A transition between phases II and 
III. In phase II dipole fluctuations of the solute relax to zero 
for long times, while they are frozen in phase III. The value 
Cc = Cc(Vc) = 0.380 is marked by an arrow. The dashed line 
is the corresponding transition curve calculated in Ref. | p9| ; 
it terminates at f' = 0.297. 



C = 0.8 




FIG. 2. Nonergodicity parameters (heavy full 

lines) for the molecule's arrested number-density fluctuations 
(x — N) and "charge" -density fluctuations (x = Z) for the 
critical packing fraction ip = ip c . The elongation parameter 
£ = 0.80 is representative for strong steric hindrance for re- 
orientational motion. The light full lines are evaluated with 
Eq. (A7) with the nonergodicity parameters f c (q£0) obtained 
in Ref. |13| from a theory based on a tensor-density descrip- 
tion. The dotted lines show the contributions to Eq. (A7) 
from different angular-momentum index I, Here and in the 
following figures the diameter of the spheres is used as unit 
of length, d — 1. 



15 



(?-C )/? = 0.347 



1 1 

1 

" 1=0 \\ 


1 1 1 1 

1 

f Nc 

1 q 

- 


1=2 ... \ 
-r--" | 






ii^-i — — ~f i 

r ZC 

1 q 


" 1=1 \ 

i 


<v 1=3 







10 20 

q 



0.5 


1 

0.5 





H — i — h 



H — i — I — i — h 




FIG. 3. Results as in Fig. ^ but for small elongations which 
are representative for weak steric hindrance for reorientational 
motion. The relative distance from the transition point be- 
tween phases II and III is (£ — Cc)/Cc = 0.347. The heavy full 
line shows the result of the present theory for £ = 0.512. The 
light full line shows the result for £ = 0.400 based on Ref. JL3J. 



FIG. 4. Critical nonergodicity parameters f% c for the so- 
lute as function of the elongation parameter £ for the wave 
numbers q = 3.4(a), 7.0(6), 10.6(c), 14.2(d), and 17.4 (e). 
The arrow marks the transition point from phase II to phase 
III at Cc = 0.380. 
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FIG. 5. Correlators ^ (i) (solid lines) and 4>q (t) (dashed 
lines) for two intermediate wave numbers q as function of 
the logarithm of time t. The decay curves at the critical 
packing fraction <p c for number-density and "charge" -density 
correlators are shown as dotted lines and marked by N and 
Z, respectively. Only a few solutions of glass states are shown 
for <j}q (t) in order to avoid overcrowding of the figure. The 
distance parameter is e = (<p — tp c )/ip c = ±10~ x . The full 
circles and squares mark the characteristic times t a and t' a , 
respectively, according to Eq. ( |36| ) for x = 1, 2, 3 and 4. The 
unit of time is chosen here and in the following figures such 
that the thermal velocity of the solvent reads v — 1. 
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log 10 t 

FIG. 6. Reorientational correlators Ce (t) for 1 = 1 (dashed 
lines) and I = 2 (full lines) for two elongations The cor- 
relators for the critical packing fraction tp = tp c are shown 
as dotted lines marked with a. The distance parameters 
e = (ip — ip c )/tp c are -0.01 (faster decay) and -0.001 (slower 
decay). The full circles and squares mark the corresponding 
time scales t a and t' a , respectively, from Eq. (|36|). The open 
circles and squares on the curves mark the characteristic time 
scales rf } and respectively, defined by Ci(rp) = f% and 

c t {?i) = ft/2. 
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FIG. 7. Susceptibility master spectra x"{&) °f the 
a-process as function of the logarithm of the rescaled fre- 
quency u) = Ljt'„ (see text). Upper panel: the curves I = 1 
and 2 refer to the response for the dipole and quadrupole, re- 
spectively, for elongation £ = 0.80. The curve I — refers to 
the susceptibility master spectrum of the dimensionless longi- 
tudinal elastic modulus m q =o(t) of the HS S. The dashed lines 
exhibit the von Schweidler tails, Eq. (39b). The dotted lines 



are fits by Kohlrausch spectra Xk(^) with stretching expo- 
nents [5 = 0.97, 0.88 and 0.63 chosen for i = 1,2 and 0, re- 
spectively, so that the maximum and the full width at the half 
maximum W in decades of x^-(ili) agree with those of x"(ui). 
The position of the susceptibility maximum is tlwx = 0.337 
(0.927, 2.69) and the width is W = 1.17 (1.28, 1.76) for £ = 1 
(2, 0). Lower panel: corresponding results for £ = 0.43. The 
stretching exponent /3, the maximum position Q max and the 
width W for I = 1 {I = 2) are /3 = 0.79 (0.71), w max = 17.5 
(2.11) and W = 1.42 (1.57), respectively. 
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FIG. 8. Double logarithmic presentation of the normalized 
fluctuation spectra of the dipole-reorientation a-processes, 
Xi (aj)a}max//ia), as function of dj/tlwx. Here £i ma x denotes 
the position of the susceptibility maximum. Following Dixon 
et al. |u| the vertical axis is rescaled by to -1 and the hor- 
izontal one by to _1 (1 + w _1 ), where w — W/Wd is the 
ratio of the logarithmic full width at half maximum W of 
the susceptibility peak to the same quantity Wn of a De- 
bye-peak. The open circles rep roduce some of the dielec- 
tric-loss results for glycerol |3l|. The three full lines from 
the top are the results for C, — 0.80, 0.60 and 0.43, succes- 
sively, although the upper two curves cannot be distinguished 
within the resolution of the figure. The dotted and the dashed 
lines exhibit the Kohlrausch fit with the stretching exponent 
/3 = 0.97 and the von-Schweidler-law tail, respectively, for the 
£ = 0.80-spectrum. 
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FIG. 9. Dipole correlators Ci(t) for the distance parameter 
e = — 10~ 4 for three elongations £. The full circle and square 
indicate the times t a and t' a , respectively, from Eq. (p6|) . The 
dashed lines are calculated from Eqs. (Gq) and (|29|). The 
plateaus /f = 0.905 (0.769, 0.376) for £ = 0.80 (0.60, 0.43) are 
shown by dashed horizontal lines. The full lines exhibit Ci(i) 
for the same states,but evaluated from equations derived 
in analogy to Eqs. (|3^)-(|35|), and their plateaus fl — 0.907 
(0.782, 0.402) are indicated by full horizontal lines. 



1.0 
0.8 
^ 0.6 



- 


? = 0.43 \\ 




\ C, = 0.80 \ 




x=2\ 


x=3\ 


x=2 \ x=3 \ 


-1.0 
0.8 






\ \ \ 


. \\ £ = 0.80 . 
\ \ x=3 
5=0.43 






0.0 t 1.0 





-4 4 8 12 




FIG. 10. Dipole correlators Ci(t) for a dumbbell of two 
fused hard spheres of diameters d and distance C,d between 
the centers moving in a liquid of hard spheres with diame- 
ter d for a distance parameter (tp — if c )/ip c = — 10~ x . The 
dashed lines reproduce the results from Fig. ^| and refer to 
a symmetric molecule with masses of the two atoms being 
equal to the mass m of the solvent particles rriA = tub = to. 
The full lines exhibit the results for an asymmetric dumbbell 
with m,A = 10m, mj = m. In the main frame the results 
for different states are successively shifted horizontally by 2 
decades in order to avoid overcrowding. The inset shows the 
transient dynamics on a linear time axis. 
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